A theoretical analysis of inertia-like switching in magnets: applications to a synthetic 

ant iferromagnet 
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The magnetization dynamics of a synthetic antiferromagnet subject to a short magnetic field pulse, 
has been studied by using a combination of first-principles and atomistic spin dynamics simulations. 
We observe switching phenomena on the time scale of tens of picoseconds, and inertia-like behavior in 
the magnetization dynamics. We explain the latter in terms of a dynamic redtstnbution of magnetic 
energy from the applied field pulse to other possible energy terms, such as the exchange interaction 
and the magnetic anisotropy, without invoking concepts such as inertia of an antiferromagnetic 
vector. We also demonstrate that such dynamics can also be observed in a ferromagnetic material 
where the incident field pulse pumps energy to the magnetic anisotropy. 

PACS numbers: 75.78.-n, 75.50.Ee, 75.50.Ss, 75.75.-c 



I. INTRODUCTION 

Studies of magnetization reversal by an external mag- 
netic field is a central paradigm which relates to mag- 
netic data storage.— The magnetic random access mem- 
ory (MRAM) device^, based on magneto-resistance (MR), 
where magnetic bits can be described in terms high (1) 
or low (0) resistance of MR units, is found to be one 
of the most promising memory devices for the near fu- 
ture. The magneto-resistance unit consists of a non- 
magnetic layer sandwiched between two magnetic metal- 
lic layers. The magneto-resistance could be either gi- 
ant magneto-resistance (GMR)2.:^ or tunnel magneto- 
resistance (TMR)^i^ depending on the sandwiched layer, 
which could be either metallic or insulating. Such mem- 
ory devices, with a very dense architecture, capable of 
performing with low power consumption is proposed to 
be a Universal memory device with advantages of other 
memory devices such as static RAM (SRAM), dynamic 
RAM (DRAM) and flash memory^. An MRAM, which is 
capable of performing at the speed of SRAM with dense 
architecture comparable to DRAM and non-volatility, 
(i.e. there is no loss of information when the power is 
switched off) is to some extent limited by it's writing 
speed. This speed depends on the magnetization rever- 
sal time, where the magnetization of a given magnetic 
unit is reversed. 

The potential for applications of magnetization dy- 
namics in technology, combined with the new knowledge 
of fundamental magnetic interactions in nano-sized ma- 
terials, has caused a large increase in the interest in this 
field. Several mechanisms for shorter switching speeds 
have been proposed, e.g. toggle switching'', all optical 
switching^ etc. One recent suggestion is the so called in- 
ertia driven switching.^ Here a short magnetic field pulse, 
optically generated, of 100 fs, was found to cause a dy- 
namic behavior of the magnetism of HoFeOs, which went 
on for tens of picoseconds. In the paper of Ref. it was 
pointed out that this behavior is expected only for mate- 
rials with antiferromagnetic exchange interactions, either 
for bulk compounds with intrinsic antiferromagnetic ex- 



change (e.g. like HoFeOa) or in artificial magnetic tri- or 
multilayers, with antiferromagnetic interlayer exchange. 
It was argued that for materials with antiferromagnetic 
exchange interactions, the switching can be described by 
a second order differential equation of the so called anti- 
ferromagnetic unit vector, and hence have features asso- 
ciated to an inertia. It was in particular argued that the 
inertia accumulated during the short pulse, was responsi- 
ble for the magnetization dynamics that lasted even after 
the pulse was switched off. 

In this paper we address the microscopic mecha- 
nisms behind this kind of switching for a magnetic tri- 
layer, using atomistic spin-dynamics simulations .-^S^li We 
have chosen for our studies a Fe/Cr/Fe model system, 
with antiferromagnetic interlayer exchange interaction 
(lEC). A schematic figure of the system considered is 
shown in FigjlJ The exchange constants were calcu- 
lated from first principles by mapping the ground state 
electronic structure to a Heisenberg Hamiltonian, us- 
ing the Liechtenstein-Katnelson-Gubanov method. We 
have calculated the time-evolution of the magnetizations 
of a Fe/Cr/Fe trilayer using spin dynamics simulation 
based on Landau-Lifshitz-Gilbert method, using the Up- 
pASD package. In our simulations we have analyzed the 
different terms that play a role in the switching behavior 
of materials with antiferromagnetic exchange, which are 
exposed to ultra short pulses. Although our simulations 
reproduce much of the observations reported in Ref.^, 
our interpretation of the microscopic mechanisms is dif- 
ferent, and we argue that a more appropriate description 
of the mechanism behind this switching is described as 
due to a redistribution of the energy from the ultra fast 
magnetic pulse to the antiferromagnetic exchange inter- 
action and magnetic anisotropy of the material. After 
the pulse is switched off the dynamics is governed by 
an effective field, now composed only of the exchange 
and anisotropy fields, that temporally evolves for much 
longer time scales than the initial pulse. The dynami- 
cally evolving magnetization carries the effect of the field 
pulse applied at the initial phase of the simulation and 
the effect persists through the course of simulation. 
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II. METHOD 

The temporal evolution of sub-lattice magnetizations 
are calculated using the dynamical equation of Landau- 
Lifshitz-Gilbert form expressed in terms of atomic 
momentsi^ 

(1) 

where mi(i) is the atomic moment on the i*'' site at time 
t. 7 is the gyro-magnetic ratio and a is the Gilbert damp- 
ing factor. The effective field H^jj- on the i*'' atom is cal- 
culated from the effective magnetic Hamiltonian given 
by 

Hmucj = Hex + HmA + Ha-, (2) 

through 

„i _ dHuag 

The sub-lattice magnetization for a given Fe or Cr sub- 
lattice is simply 

M(t) = ^mi(i). 

i 

The first term in the Hamiltonian describes the inter- 
atomic exchange interaction between between all atoms 
of the simulation cell. It can be separated into two parts; 
(i) pair-exchange between the atomic moments in the 
same layer (intra-layer exchange coupling) and (ii) the 
exchange coupling between the moments residing in dif- 
ferent layers (inter-layer exchange coupling). The second 
term of Eqn.2 represents the magnetic anisotropy. The 
form of magnetic anisotropy Hamiltonian we have con- 
sidered is Hma = Ksin^0, where 9 is the angle which 
the sublattice magnetization makes with anisotropy axis. 
K is the anisotropy energy density and can be written 
as the sum of two contributions K — Kmca + K shape- 
Furthermore, Kmca is the energy density which corre- 
sponds to magneto-crystalline anisotropy and K shape cor- 
responds to the energy density for the shape anisotropy 
arising from dipole-dipole interaction. We choose K such 
a way that the hard axis is considered to be along the 
z-direction and the xy-plane is considered as easy plane, 
due to the shape anisotropy of thin magnetic films and 
due to the possible intrinsic uniaxial anisotropy. The last 
term of Eqn.2 is interaction due to the applied magnetic 
field. For the present case it describes the pulsed field we 
apply during a short time of the simulations. Equation 
1 was solved numerically using the procedure of the Ref. 

M 

For T= K we consider the macrospin approximation, 
i.e. all atomistic spins within one layer are assumed to 
co-rotate. For finite temperature simulations, one has to 
deal with this problem in an atomistic way, considering 



a large number of atoms in the simulations^ii As men- 
tioned, the exchange constants were calculated from first 
principles by mapping the ground state electronic struc- 
ture to a Heisenberg Hamiltonian. The ground state elec- 
tronic structure was calculated using the Korringa-Kohn- 
Rostocker Green's function method within the atomic 
sphere approximation . ^^'-'^^ In reality the interlayer ex- 
change between the Fe layers through the Cr spacer layer 
are influenced by the interface roughness of the Fe-Cr in- 
terface, and to simulate this we rescaled our interlayer 
exchange interactions with a factor of 10, which is ac- 
cording to previous analysis. ^^'^^ 

We consider each Fe sublattice has magnetization M = 
Msn(t), where h{t) is the unit vector. The magnitude of 
the effective anisotropy field is given by Fma = ^cosO. 
Ms is the saturation moment for the Fe sublattice. The 
initial value of a Fe sublattice magnetization is along 
the positive or negative x-direction, as shown in Figll] 
i.e. M = ztMsBx. For all simulations we used a square 
shaped pulse of the external field. At time t=0, we apply 
a field of constant value in the x-z plane with constant 
value described by Ha = [H^ H^]. We studied the cases 
where the magnetization was in-plane, due to the shape 
anisotropy, without any further in-plane anisotropy. 




FIG. 1: (Color online) (Color online) The considered geom- 
etry of the atomistic simulations. The red slabs represent 
Fe and the blue slab Cr. The arrows in the slabs represent 
directions of magnetization. 



III. RESULTS 

We start by showing in Figl2]an example of the magne- 
tization dynamics of the tri-layer in FigHJ Note that the 
magnetization shown in the figure is only for the top Fe- 
slab in FigH] (normalized to unity). The results shown 
are for a field pulse with Ha;= 3 T and H^^S T, for a 
duration of 5 picoseconds (ps). The Gilbert damping 
constant used was 0.02 and the strength of the uniax- 
ial anisotropy constant was 0.06 mRy/atom, providing 
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FIG. 2: (Color online) The time-evolution of the x-, y- and 
z-component of the top-most Fe sublattice magnetization of 
the tri-layer shown in Fig[T] We have used }ix=3T, H2=5T, 
|K|=0.06 mRy for a = 0.02 



an easy plane magnetization. At time t=0, we applied 
a field pulse. From the figure, we see that as soon as 
the pulse is applied, the magnetization is brought out 
of plane, and all the components of magnetizations are 
non-zero. The magnetization then precess under the in- 
fluence of the external field, as well as the internal fields, 
provided by the anisotropy field and the exchange field 
due to the lEC. From the figure it is clear that even af- 
ter the external magnetic field is switched of (i.e. after 
5 ps) the system continues to evolve in time, and the 
magnetization dynamics continues over the entire time- 
interval shown in the figure. At sufficiently long simula- 
tion time it seems however that the x-componcnt of the 
magnetism approaches the negative x-direction. Hence 
for this choice of external pulse a switching of the mag- 
netism occurs, and the direction of all magnetic moments 
illustrated in Fig[T] have become reversed. The data in 
Fig. [5] are similar to those demonstrated experimentally 
in Ref. [^, in that the magnetization dynamics continue 
long after the external pulse has been first applied and 
then removed from the system. 

By sweeping the strength of the x-component of the ap- 
plied field we can realize different situations, in which the 
sublattice magnetization of one of the Fe slabs in FigU] 
rotates half a turn (180 degrees), a whole turn (360 de- 
grees) or even more than one turn. In Fig|3]we illustrate 
this result, in a phase diagram of the switching behaviors. 
In the simulations shown in Fig 13] we have used the same 
z-component of the applied field, the same Gilbert damp- 
ing factor, pulsing time, and the same uniaxial anisotropy 
constant, as used in Fig (21 The figure only contains infor- 
mation about the Fe magnetization of the top Fe slab of 
Fig IT] In the phase diagram of Figl3]the time of the sim- 
ulations is shown on the y-axis and the magnitude of the 
x-component of the applied field is shown on the x-axis. 
The switched and unswitched regions are represented in 



terms of a color coding. The red color specifies that the 
magnetization has turned n ■ 360 degrees (n=0,l,2,...), 
whereas blue color specifies that the magnetization has 
turned 360 -(-180 degrees (n=0,l,2,...). Taking as an ex- 
ample Hx— 5T, Figl3^ shows that when the pulse arrives 
the magnetism of the top Fe slab is along the positive x- 
direction (red color), but after 5-6 ps the magnetization 
has reversed to the negative x-direction, and after 10 ps 
a rotation of 360 degrees to the positive x-direction has 
occurred. In between, the magnetism of course points 
along the positive or negative y-direction (this compo- 
nent is not shown), and hence exhibits a rotational be- 
havior. After some 30 ps the magnetism finally settles 
down along the negative x-direction. Overall, the data 
in FigO shows that the magnetization dynamics contin- 
ues long after the external pulse has been turned off, in 
a similar fashion as observed experimentally.— 

Another interesting feature of the data in Figl3^ is that 
the results are not symmetric with respect to positive or 
negative x-component to the applied field. This might 
at first sight seem puzzling, since the geometry shown 
in FiglT] is very symmetric with respect to positive and 
negative x-direction. However, the data in Figl5^ depict 
only the magnetization of the topmost Fe layer, and for 
this layer alone there is no obvious symmetry for cases 
when the applied pulse is along the positive or negative 
x-direction, and the patterns in Figj3^ for positive and 
negative values of H^; are different. However, for all val- 
ues of Hj; we find that after sufficiently long time the 
magnetism has settled down into a given orientation. 

There are two torques that in general are relevant for 
the magnetization dynamics, a precessional torque and 
a damping torque, and the influence of the latter is pro- 
vided in the simulations by the strength of the Gilbert 
damping parameter. Hence, one may expect that the 
phase diagram should depend critically on this parame- 
ter and to investigate this possibility, we show in Fig. JJb 
a simulation with all parameters kept exactly the same 
as for the simulation shown in Figj3^, with the one dif- 
ference that the Gilbert damping parameter was set to 
0.002, i.e. ten times lower than that used for FiglS^. 
In this case the rotation of the magnetism continues 
throughout the entire period of the simulations, at least 
for positive values of H^; , which shows that it is the damp- 
ing parameter that primarily determines when the mag- 
netization stops rotating. 

As the pulse hits the sample (simulation box) the an- 
gle between the two Fe slabs start to deviate from the 
180 degree antiferromagnetic coupling. To illustrate this 
fact we show in Fig|3] a contour plot for the angle be- 
tween two Fe sublattice magnetizations, in which red col- 
oration signifies 180 degree antiferromagnetic coupling. 
This plot corresponds to the simulations shown in FigI3^ 
and Figl5}D. As is clear from Fig|3^,b, the magnetiza- 
tion of the two Fe sublattices are initially at 180 degree 
with each other, but after the external field is switched 
on and time evolves, the angle between the sublattices 
is reduced from 180 degrees and oscillates with time, es- 
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(b) For a=0.002 
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FIG. 3: (Color online) The switching phase diagram of an 
Fe/Cr/Fe-trilayer at T=0 K. The duration of the magnetic 
pulse is 5 ps, and the x-component of the magnetic field is 
shown on the x-axis, whereas the z-component is 5 T. A uni- 
axial anisotropy of 0.06 mRy is considered in the calculations, 
providing an easy plane magnetization. The x-component of 
the magnetization of the top Fe slab of Fig[l] is shown in a 
color scheme. The blue/red areas indicate the regions where 
switching has occurred/not occurred (see also text for de- 
scription). In a) we show simulations for a damping constant 
of 0.02 and in b) the damping constant is 0.002. 11^= 5T, 
|K|=0.06 mRy 



pecially when H^; ^ 0. Fig|3^,b shows data for two dif- 
ferent damping parameters, and we note that these two 
different damping parameters result in slightly different 
magnetization dynamics, in that it takes longer time for 
the angle between the two sublattices to approach 180 
degrees for the weakly damped system. However, for 
both simulations a 180 degree coupling is obtained when 
t — >■ cx). The minimum angle is seen when the magnitude 
of the x-component of the applied field is largest, after 5 
ps (i.e. the duration of the pulse), where the two sublat- 
tice magnetizations are approximately perpendicular to 
each other. This means that, with the applied field as a 
driving force, the angle between the two sublattices is in- 
creasingly deviating, albeit in an oscillating fashion, from 
antiferromagnetic coupling. When the angle between two 
sublattices is less than 180 degrees an internal exchange 
field is build up in each Fe slab, which causes a dynam- 
ical response. In addition, a non-zero z-component to 
the magnetization (e.g. as shown in FigI2]) builds up an 
anisotropy field. These two contributions to the effective 
field, which drives the magnetization dynamics, are non- 
zero long after the external field has been switched of. In 
fact, they are zero only after the magnetization has set- 
tled down to a constant value. Hence, the effective field 
is in this case dynamic, and the torques due to this field 



has as usual two components; the damping part and the 
precessional term. 

The data in Fig]?] demonstrate an important result, 
namely that during the short time that the external field 
is applied, a significant energy is established both in 
the exchange interaction term, i/ea;, as well as in the 
anisotropy energy, Hmae, of Eqn.2. In a way one can 
view the whole process as being due to energy provided 
by the Ha-term of Eqn.2, that during the pulse duration 
pumps energy into the other two terms, Hg^ and Hmae- 
Even if Ha is non-zero only for a short time, the other 
two terms will be non-zero for a much longer period, and 
the effective field composed of contributions due to these 
two terms depending on mi(t), carries some sort of in- 
formation of the field pulse applied at an earlier stage. 

To demonstrate this in a more formal way, we plot in 
Figl5]the temporal evolution of the different components 
of magnetic energy during one of the simulations. From 

(a) For a=0.02 
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FIG. 4: (Color online) Evolution of angle between the sublat- 
tice magnetizations (a) for damping factor of 0.02 (b) with a 
damping factor of 0.002. 5 T, |K|=0.06 mRy,a=0.02 

FiglS] one can see that the Zeeman energy is negative 
during the duration of the pulse whereas the exchange 
and anisotropy energy are positive. The sum of the three 
components, i.e. the total energy is negative during the 
duration of the pulse, which means that during the time 
that the pulse is applied, the energy is lowered and the 
system tries to relax into a new ground state configu- 
ration. When the external pulse is switched off, FiglS] 
shows that there is still positive energy in the exchange 
and anisotropy contribution to the energy and the total 
energy is now positive. In this configuration, there is both 
a z-component of the magnetization as well non-collinear 
coupling between the moments in the two Fe-slabs. The 
system then relaxes back to the original configuration of 
coUinear antiferromagnetic coupling of the Fe-slabs with- 
out a z-component to the magnetism and it does this in 
the anisotropy and exchange field that has been build up 
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during the duration of the external field . The relaxation 
back to collinear configuration takes considerably longer 
time than the duration of the field as FiglS] shows. From 
FiglS]it can be seen that the energy due to applied field 
(Zccman energy) acts as a reservoir for the other two 
terms (exchange and anisotropy) of the magnetic Hamil- 
tonian. Hence, in order to analyze the magnetization 
dynamics provided in FigsH] and Fig. [31 one does not 
need to invoke concepts like inertia of antiferromagnetic 
unit vectors^. We suggest a simpler interpretation based 
on a redistribution of energy among the three terms in 
Eqn.2, in particular a redistribution of energy from Ha to 
Hex and Hmae , during the time period that the external 
field is applied. 
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FIG. 5: (Color online) Time evolution of different components 
of magnetic energy (Exchange, Anisotropy, Applied), plotted 
alongside with total magnetic energy. Hz= 5 T, a=0.02 



we note that also in this case can the magnetization be- 
have as if it had 'inertia', i.e. it proceeds after an external 
field is switched off. An analysis relying on antiferromag- 
netic unit vectors becomes cumbersome in this case, but 
the pumping of energy from Ha to the two other compo- 
nents given in Eqn.2, in this case the anisotropy energy, 
provides an easy interpretation. The data in Fig. 7, are 
actually consistent with the experiments performed in 
Ref. Il9l . where a short pulse was applied to a ferromag- 
net thin film. 




0.02 0.04 0.06 0.08 0.1 
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FIG. 6: (Color online) Switching times as function of out-of- 
plane anisotropy. The duration of the magnetic pulse is 5 ps, 
the z-component of the applied magnetic field is 5T and the 
x-component is 3T. (All other parameters are same as Fig.2) 



The analysis in the preceding paragraph indicates that 
also the strength of the anisotropy influenced the mag- 
netization dynamics and to investigate this further, we 
plot in FigEl simulated data as a function of uniaxial 
anisotropy constant, K. All other parameters were the 
same as for the simulation of Figl3^. First we note that 
when K=0, the tendency of switching is reduced, since in 
this case the dynamics is provided by the exchange field 
alone, via the lEC. However, when K increases to 0.01- 

0. 04 mRy, the tendency for switching increases and the 
magnetization is seen to precess two times before stabiliz- 
ing along the positive x-direction. For even larger values 
of anisotropy energy, K larger than « 0.05 mRy, only one 
revolution is found in the simulations, before the magne- 
tization is stabilized. Hence there is a critical value of K 
which provides an optimally large amount of revolutions 
in the magnetization. 

We have also investigated the case when the inter- 
layer exchange coupling is zero, in Eqn.2, by perform- 
ing simulations to only the top-most Fe slab, of Fig. [1] 

1. e. a purely ferromagnetic system having an easy-plane 
anisotropy. The simulated results are shown in Figl7]and 
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FIG. 7: (Color online) Switching phase diagram for Fes lay- 
ers with the same set of parameters as indicated in Fig. 3 for 
Fe3/Cr4/Fe3 trilayer, except the lEC. Hz= 5 T, [K|=0.06 
mRy,a=0.02 
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IV. CONCLUSIONS 

In conclusion we have performed atomistic spin dy- 
namics simulations for a synthetic antiferromagnet sub- 
jected to a short magnetic pulse. We have demonstrated 
that the magnetization dynamics continues long after the 
pulse has been applied and removed, in a fashion that is 
consistent with the experimental results of Ref. 0. 

Several materials parameters, like the damping and 
strength of the magnetic anisotropy, are found to in- 
fluence the details of the magnetization dynamics. We 
suggest that such magnetization dynamics could be un- 
derstood in terms of a redistribution of energy, during 
the time that the pulse is applied, so that the magnetic 
anisotropy energy and interlayer exchange energy become 
in an out-of-equilibrium state. The magnetic energy due 
to the incident pulse hence becomes distributed dynam- 



ically to the effective interlayer exchange coupling be- 
tween the Fe layers and it also modifies the magnetic 
anisotropy energy. The latter term ensures that the ob- 
served behavior of a magnetization dynamics that pro- 
gresses long after an external field has been applied and 
switched off, is possible also in ferromagnetic materials. 
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